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f I : Abstract. Several recent experiments have demonstrated the viability of a passive device 

f~*i ■ that can generate spin-entangled currents in two separate leads. However, manipulation and 

^^ , measurement of individual flying qubits in a solid state system has yet to be achieved. This 

is particularly difficult when a macroscopic number of these indistinguishable qubits are 
present. In order to access such an entangled current resource, we therefore show how to use 
k> \ it to generate distributed, static entanglement. The spatial separation between the entangled 

Vh ' static pair can be much higher than that achieved by only exploiting the tunnelling effects 



between quantum dots. Our device is completely passive, and requires only weak Coulomb 
interactions between static and flying spins. We show that the entanglement generated is robust 
to decoherence for large enough currents. 
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1. Introduction 

Entanglement is an enabling resource for quantum computing (QC). It must be created and 
consumed in the process of executing any quantum algorithm fl}, something which is most 
obviously apparent in the measurement-based model of quantum computing [|2l[3]|. In this 
picture, entanglement is first generated to build cluster states, before being consumed by 
single qubit measurements during the execution of an algorithm. The initial entanglement can 
be generated between distant qubit nodes HIEIO, each of which can have its own, dedicated, 
measurement apparatus. Distributed entanglement would also enable secure communication 
over long distances D7l[8]| and quantum teleportation (HJ. 

Devices that generate entangled currents of pairs of electron spins propagating 
down different leads have been proposed in theoretical work IfTOl ITTl [T2| . and recent 
experiments lfT3l [T4ll have begun to demonstrate their feasibility lfT51 [T6l [TTl . However, 
it is not clear how such an entangled resource could be used for any of the applications 
discussed above, since the control and measurement of a single flying solid state qubit 
has yet to be demonstrated experimentally. Furthermore, the macroscopic nature of the 
currents makes this even more difficult, especially when there are times when spin pairs 
enter the same lead |[T3l[T4l . In this paper, we will show that it is possible to convert such 
mobile entanglement to a static form in a completely passive way, in a very simple device 
- thus opening up the possibility of quantum information processors based on entangled 
currents. Notably, our scheme produces a static pair of entangled electron spins that can 
have a much higher degree of separation (see Sec. 14.21) than more conventional protocols for 
entanglement generation for which the separation is limited by quantum tunnelling or similar 
local interactions lfT5lfT6l[T7ll . 

At the centre of our device is a spin-entangled current generator that outputs entangled 
pairs of spins fa down two different leads i. Each spin encounters a further, static, spin s, 
downstream of the generator and interacts with it, as shown in Fig. [Q The generator is based 
on earlier proposals of a passive device that produces pairs of spin-entangled electrons, with 
each pair maximally entangled in the singlet Bell state \S)f = 4=(| t/il/ 2 ) ~~ I 4/it/ 2 )) 
IfTOl [TTl [T2l . The nature of the static spins is not important, but one possibility is that a 
single electron is confined in each of the two quantum dots that are fixed close to the leads. 
Other suitable architectures include endohedral fullerenes in carbon nanotube peapods |[T8l . 
carbon nanobud structures |fT9ll , and surface acoustic waves whose minima isolate single 
electrons l|20l . 

2. Model and Basic Results 

Let us start with an effective Hamiltonian coupling the flying and static spins of the following 
form: 

Hi = f {o"c{ + ofcr}) = 9i (a s i<r f j + <£erj) (1) 

where the a± = (a x ±ia y )/2 are the usual Pauli operators. The ^ are XY exchange coupling 
strengths that depend on the separation of the two spins. Each g^ is time dependent since one 
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Figure 1. Illustration of our entanglement generation device. Two currents, composed of 
successive electron pairs that are maximally spin-entangled, emerge from a generator and pass 
down two leads. A static electron is situated downstream of the generator close to each lead, 
and the static pair is spatially well separated. The mobile spins in each lead interact with 
corresponding static spins as they pass. Note that the two flying spins of the same Bell pair are 
not required to, and will normally not, arrive at the sites of interaction with their corresponding 
static spins at the same time. 



of the two interacting spins is mobile. The time evolution operator Ui(t) for a general state 
\^i(t)) = Ui(t)\^i(0)) of the static-mobile pair i, is then 

Ui(t) = exp [-^(£)(ctJct^ + criv^)] , (2) 

where 9i(i) = f Q gi(t')dt' /h. 9i(t) is constant when [0, t] is chosen so that ^(O) and gi(t) are 
negligible. In the basis | t*l7«). I UW> I Ut/*), I UW> 

( 1 \ 

cos 9, i —i sin 9. h 

—i sin 9i cos 9, t 

1 



Ut 



\ 



(3) 



/ 



2.1. Attractor 



We first consider the case where 9\ — 9i — 9 and initially the static spins are in the 
state | tsits 2 )- The starting state of the two static and the first pair of flying spins is then 
TTiO TsiT/i)l Ts 2 4-f 2 ) — I tsil/i)| Ts 2 Tf 2 ))- Using the product of 2 unitary operators of the form 
of Eq. [31 we find that after interaction the state becomes cos 9\ X\) s \S)f — « sin f5|S'} s | ft) f- 
Since this first flying pair will no longer interact with the static spins nor with the following 
flying pairs, we can trace out this first flying pair to find the density matrix describing the 
static pair after the interaction: cos 2 9\ TT)s(TT | + sin 2 9\S) S (S\. 

In order to find the behaviour of our system for multiple passages of flying qubits, we 
require the following four maps which describe a single interaction event, C/i, Ui'. 

I TTMTT I H- cos 2 9\ tt) s (tt I +sin 2 f9 \S) S (S\; 
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I U)s(Ur I >->• cos 2 9 I U) s (ll \+sm 2 9\S) s {S\; 

\S) S (S\ M- ^sin 2 2#(| tt) s (tt I + I ii)s(ii \)+sin 4 9 |T ) S (T | + cos 4 # |S) S (S|; 



|To) fl (To| ^|To) s (T |, (4) 

where |T ) = 4=(| tl) + I It})' another maximally entangled Bell state. These maps imply 
that the state |T ) S is an attractor for this process, and by tracking the states of the static spins 
through multiple passages of the flying spins it is clear that the system will converge towards 
this maximally entangled attractor; we need not consider any further maps since other static 
spin states are never accessed. This fixed point is also consistent with spin-invariant scattering 
of the flying singlets from a static \T ) triplet, as can be verified directly from the Schrodinger 
equation. 



2.2. Convergence 

We now calculate the probabilities of obtaining the \T ) S state after a number of flying spin 
passages. The reduced density operator for the static spins after n iterations is, by definition, 



» 



is the vector 



pi' = P re ■ &, where & = (J tt),(tt I, I U)s{U I, \S) S (S\, \T } S (T 

of projection operators for the base states and P„ = (P^(n), -P^(n), Psiji), r T 
the corresponding vector of probabilities. Under the map, Eq. & H- L3? and hence 

J»,i.e., 



rp 

PrJn)) is 



>) ^ > +1 ) 

Ps ^r Ps 



p n . (l <?) = (L T p n : 

P„ ^ P n+1 = M P„ 
where, by direct substitution, 

/ cos 2 e 
o 



(5) 
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sin 2 6> 
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sin 2 9 





\ sin 2 29 

\ sin 2 29 

cos 4 # 
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Note that the map, Eq. |51 preserves total probability, as it should. Since the matrix M 
describes the evolution of the state probabilities, its eigenvalues A« must satisfy |A»| ^ 1. 
Except when 9 is a multiple of n, there is always one (and only one) eigenvalue equal to unity, 
and the corresponding eigenvector is then our attractor state |T )(T |. Using Mathematica 
and Eq. [51 we have derived closed analytic expressions for M n and P To [n). For the initial 
conditions P = (1, 0, 0, 0) T , we have 

5 + cos 29 



PT (n) 



-ira-l 



[{B n + C n ) 



A 



■(B n -C n )} 



(6) 



where 



.4 
B 
C 



V2 cos 2 9(17 + cos 29) 

7 + 8 cos 29 + cos 49 - AA sin 2 9 

7 + 8 cos 29 + cos 49 + AA sin 2 9. 



(7) 
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For large n, when (4Asin 2 6)/ (7 + 8 cos 29 + cos 49) e (0,2), i.e., when 9 < 0.89, C 
dominates over B n and we have 



ln(l-P To (n)) ~a(9)n + (3(9) 



(8) 



where a{9) = ln(C/16) ^ and (3(9) = \n[(A + 5 + cos29)/(2A)}. For Q < 1, 
(5 + cos29)/A -> 1, and (3(9) ->■ 0. We plot the probability of obtaining the |T ) S (T | 
state against 9, for different values of n in Fig. |2(a)[ 
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Figure 2. Plots of the probabilites of obtaining the \Tq) s (Tq\ state (a) against -, after 
10, 50, 100, 500 rounds of interactions, respectively; (b) - (d) against the number of rounds 
of interactions, n, for 8 — 0.1, 0.03 and 0.01, respectively. 



As the number of rounds increases, the range of 9 values for which there is a high 
probability of obtaining the static \T ) state becomes wider; we plot the probability Pr (n) 
against the number of rounds, n, for some weak coupling strengths in Fig. [2£b)-(d). For the 
tunnelling rates T in |fT4]|. the time interval t between successive flying qubits are on the order 
of f — 67 ps, and therefore the time it takes for the static spins to converge to the \T ) state 
would be 20 /is, 2 ms and 0.2 s for 9 = 0.1, 0.03 and 0.01, respectively. These time intervals 
to are at least one order less in lfT3l . and have the potential of being shortened further. Within 
the electron spin coherence time (> 200 fj,s) observed in I12TI . the convergence can occur for 
9 as small as 0.03 with the tunnelling rates achieved in lfT3l . We also point out that molecular 
systems have the potential for phase shifts 9 > 1 due to larger exchange interactions arising 
from nanometer length scales; for example the exchange coupling in a nanotube/fullerene 
system may be several orders of magnitudes greater than in gated semiconductor devices |[2"2~1|. 



Generating distributed entanglement from electron currents 6 

3. Generalisations 

Let us now generalize our analysis, to include arbitrary coupling strengths (9\ ^ 9 2 ), and 
arbitrary starting states for the static qubits. With the time evolution operators £/j defined as 
in Eq. [3l we can find a completely positive map C s which represents the effect on the static 
spin density operator p s of a passing flying qubit pair: 

p^ = C.[pP] = Tr f [(U 1 (g)U 2 )(pW(g)pf ) )(t/i(g)C/ 2 ) t ] (9) 

where Try denotes the partial trace [QQ| over the mobile spins, and pf = \S)f(S\ . 

Eq. [9] corresponds to a set of 16 recurrence relations for the elements of pi . When 
Oi = 9 2 , four of these relations decouple from all the others, as we found in our 
argument earlier. The superoperator C s is not a linear map, and we thus vectorize the 
density operator states by listing the entries in the following order as a column vector: 

(pnPi2---Pi4P2i---P24P3i Paa) T ='■ P- The map corresponding to the action of the 

superoperator C s , 

L : p{ k) H4 p? +1) (10) 

is then linear and can be written as a simple 16 x 16 matrix M, whose entries can be easily 
calculated using Eq. |9j 

We find that M always has an eigenvalue A = 1 , independent of the values of the coupling 
strengths. The multiplicity of this eigenstate is one, unless 9 1 and 9 2 are multiples of n. The 
corresponding eigenvector pi is then a single attractor state that is independent of the initial 
configuration of the static spins, which when transformed back to its density matrix form is 

/ a \ 



fh 



2(a + b) 
where 



h 75 ° 
71 b ° 



(11) 



\ a ) 

cos #1 — cos # 2 ) 2 (1 + cos Q\ cos 8 2 ) esc 3 Q\ esc 3 9 2 
2^2 



>0, 



esc 6 1 ! esc ^2 (esc 2 ^1 + esc 2 #2) — cot 6 1 ! cot ^(2 + CSC 2 ^! + esc 2 # 2 ) 1 

27! -72' ( } 

When Q\ = 9 2 , we have a = and b = 4=, and p\ reduces the \T ) (T \ state as in the simple 
case. When 9 1 , 9 2 E (0, 7r/2], all the other eigenvalues are numerically in the range (— 1, 1), 
and their corresponding eigenvectors, when back in matrix form, all have trace zero and hence 
do not correspond to physical density operators. This has to be the case, as can be seen from 
the following argument. We can express any density state as p = J2i=i a ipii since the set of 
eigenvectors p, (vectorized pj form a basis for M. The coefficients a { can take any values, so 
long as Tr(p) = 1 and p is positive semi-definite and Hermitian. Hence, C s [p\ = J2lt± a>iKPi, 
or more generally £™[p] = X^=i a iKPi> which also has trace 1 as a density matrix. As 
n — > 00, A" — > Vj Ai| < 1 which is when i — 2, 3, ..., 16, and thus £™[p] — > a\p\. So, 
Tr(£"[p]) — ¥ a\ since we defined Tr(pi) = 1, and this requires a x = 1. We thus obtain 
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p = pi + S^£ 2 a jP«' an d the trace requirement result in £*£ 2 ajTr(pj) = 0, which holds for 
various combinations of o^'s. This can only be true when Tr(pj) = Vz = 2, 3, ..., 16. □. 
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Figure 3. Contour plots of (a) the fidelity between the attractor p\ and the |To)(To| state; 
and (b) the entanglement of formation Ep of p\, each are plotted as a function of 0\ and 62 
G(0,7r/2]. 



Now, we can find how close p\ is to the |T ) (T 1 state, for various values of 9\ and 9 2 , by 
calculating the fidelity F as defined in ll2~3~1 . In our case, we have 



F(p u \To)(T \) = Tr Jy/pi\T )(T \Jpi 



1 



(13) 



2(a + b) 

which takes a value of unity when 9\ = 9 2 , as expected. Its contour plot in Fig. |3(a)| illustrates 
that even when 9\ and 9 2 are different, p 1 is still very close to the |T )(T | state for any 
(9±, 9 2 ) in the central region. The fidelity values also indicate the levels of degradation in our 
entangled resource through unequal coupling, the degree of which can be further established 
by calculating the Entanglement of Formation IJ24l E F that the bipartite state p\ has [25]. We 
construct the contour plot for Ep in Fig. |3(b)| showing that the degree of entanglement p\ 
possesses is very large, higher than 0.9, for any (0 l5 9 2 ) in the central red region. 

The attractor state p\ remains similar for (9i,9 2 ) values in the same colour region 
in Fig. [3l and close to |T )(T | for the central red region. As a result, a high degree 
of entanglement with E F > 0.9 is achieved for the static spins even when 9i(k) varies 
with the iteration k, as long as in each round A9 = \9\ — 9 2 \ satisfy the condition set in 
Fig. |3(b)| such that (6^, 9 2 ) lies in the central red region. The maximum possible percentage 
difference in coupling strengths needed to maintain a highly entangled state may be defined as 
k(9i) := max EF > a9 (A6>/f9 1 ), and this quantity depends on 9 { itself, as shown in Fig. |4(a)[ The 
maximum value for k{9i) occurs at a #1 value away from n/2 because the maximum allowed 
A9 values for final Ep(p s ) > 0.9 do not vary much for 9{ close to ir/2. For 9% as weak as 0.1, 
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Figure 4. (a) Plot of the maximum allowed percentage difference k{Q{) against 9i/ir for 
obtaining final Ep(p s ) > 0.9; (b) & (c) Example plots of the evolutions of Ep for the static 
spins as a function of the number of flying pairs n, for strong and weak random couplings, 
respectively. The inset in (c) zooms in on the corresponding fluctuations. 



the maximum allowed percentage difference k(0.1) — 1.6%. We show two example plots for 
the evolutions of Ep for the static spins vs n for 9i values that can vary within some range 
from round to round, in Fig. |4(b)| -(c). From these results we see that the convergence times 
are of the same order as for the case when 9 1 = 9 2 (see Fig. 2) for both strong and weak 
couplings. We can also see that the condition on k{9i) in each round as specified in Fig. |4(a) 
only needs to be satisfied by the majority of the interaction rounds. 



4. Decoherence 

In practice, the environment interacts weakly with both static and flying spins, causing their 
decoherence. This will degrade our final entangled resource to some extent, and we shall 
consider such effects in this section. 



4.1. Static 

We first take into account the decoherence of the static spin pair, which in general will be 
coupled to the surrounding environment. The process can be modelled by a Markovian master 
equation |f26l with decoherence channels k and corresponding decay rates 7*.: 



,15 

J+P^t) = ^lk(A kPs A\ - -A\A kPs 



fc=i 



\psA\A k ) 



(14) 



where we have ignored the unitary term —i[H, p s ] on the assumption that the duration of the 
spin-spin interaction is much shorter than time interval t between successive flying qubits. 

This assumption about the much shorter timescale also means that we do not need to 
consider decoherence during interactions events. Therefore, the evolution of p s (t) can be 
tracked through application of the map L s to describe interaction events and using the solution 
of the master equation between events. The resulting behaviour is a function of the products 
"fkto, and we assume to to be constant for simplicity. 

For certain decoherence models, analytical results can be found. For example, if we 
assume that there are two independent dephasing channels erf 1 (£) P 2 and I Sl erf 2 with rates 
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7x and 72 respectively, then for 9\ = 9 2 the final fidelity of the state of the static spins with 
respect to |T ) is, immediately prior to an interaction events: 



2 + 2e 2 (^+^)*o csc 4 _ cot 2 0(2 + 2 esc 2 9) 
4 e 2(7i+7 2 )io esc 2 9(2 esc 2 6-1) - 4 cot 2 5(1 + 2 esc 2 



(15) 
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Figure 5. Plots of the state of the static spin pair in the presence of pure dephasing errors 
at random 0* € [1.3, 1.4] and 71 = 72 = 1 MHz: (a) Fidelity of final states at equilibrium 
against current with the higher and lower lines representing the states right after and before 
XY couplings, respectively; (b) Ep against time for a spin-entangled current of 10 pA, with 
the dashed and solid lines representing the initial static spins being in the \Tq) and the pure 
I 7t) states, respectively. For the latter case, the initial rise and the subsequent fall in Ep is 
due to the presence of the singlet state. 



We know that for the device to generate a large degree of static entanglement, we need 
Q\ and 9 2 to be in the central red region of Fig. |3] For example, working with random 
9i G [1.3, 1.4], calculations show that robustness to decoherence requires that 7 fc t are at 
most on the order of 1CT 2 . For 71 = 72 = 1 MHz (i.e., a coherence time of 1 /xs), Fig. |5£a) 
illustrates the behaviour of the fidelity against current of flying pairs f-, with e the electronic 
charge. When the current is large enough, more than 10 pA (i.e., t < 10 ns) in this case, 
the fidelity only fluctuates slightly. Fig. Ob) then shows that the final state of the static spin 
pair possesses a large degree of entanglement and does not depend on the initial set up. A 
similar effect is observed for the bit-flip errors, with slightly reduced final fidelity. This 10 
pA is much smaller than the nA scales observed in lfT3l[T4l . i.e. the strong coupling regime is 
robust to decoherence of the static spins. For weak couplings, a higher current is required for 
the same decoherence rates, e.g., at ~ 100 nA (t ~ 1 ps) is needed for 6>j ~ 0.03 (or at ~ 1 
nA with spin coherence times in ll2~l~T0 . We see that decoherence will reduce the fidelity and 
hence the degree of entanglement of the steady state for the static spins, and thus in Fig.[3]the 
regions for high fidelity and high degree of entanglement will be narrower in the presence of 
decoherence. 
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4.2. Flying 

Next, let us consider the decoherence of the flying singlets before they interact with the static 
spins. The whole analysis in Sec. [3] applies when we refine our model by allowing random 
error deviations e^ of p\ from the singlet state in each round such that 

P f ] = (l- e ^)\S) f (S\ + J2^4\ S )f( S \4^ ( 16 ) 

where the errors e^ = 2~2j=i 4 f° r eacri & are due to spin-orbit coupling and interactions 
with the environment as the flying qubit propagates. The cr/s are the Pauli o"i 2 ,3 = &x, y ,z 
matrices that correspond to different errors on one of the two flying qubits. We find that 
the single attractor state now has small error-dependent terms c for the corner entries on the 
minor diagonal of p\ in Eq. [TH and both a and b now also depend on e-- '. As a result, the 
final Ep(pi) is reduced, and the larger the errors the smaller it will be. We also find that 
the strong coupling regime can tolerate larger errors compared with weak coupling. This is 
because the static entanglement built up per round of interaction is much smaller for weak 

(k) 

couplings and the errors e:- can drastically reduce the accumulated entanglement. In either 
case, the average of e^ mainly determines the final E F for the static spins for fixed 6i. For 
strong couplings, the error tolerance on e^ is of order 0.01, which corresponds to flying pairs 
with an average E F of ~ 0.8 — 0.9. This results in a travelling distance on the millimetre scale 
for the flying spins at the Fermi velocity with typical coherence time of microseconds, before 
their interactions with the static pair. This error tolerance is much smaller for weak couplings, 
but this could also be feasibile experimentally since couplings become weaker at high carrier 
speed, which also means a smaller interaction of the flying qubits with the environment before 
they arrive at the static qubit sites. 

Therefore, for a given separation between the static spins, e.g., 1 mm (3> 1 pm, of 
conventional separations [Tf5l[T6l[T7l ). the time it takes the flying electrons (travelling at ~ 10 5 
ms _1 ) to arrive at the static qubit sites is 10 ns, much less than a typical spin coherence time 
(/is). Separations of centimetre scales can be achieved for longer spin coherence times EQ. 

5. Splitting Efficiency 

We now consider the efficiency of successful splittings of the Cooper pairs via the double dots 
in the generator |[T3l[T4ll in Fig. 1. When unsuccessful, a singlet pair enters the same lead and 



this has the effect of reducing the static entanglement (see Fig. |6(a)[ ). For a success rate p of 
50% in |[T4| . the static spins converge to a completely mixed state regardless of the coupling 
strengths. We therefore classify the usefulness of the entangled current resource depending on 
the success rates in terms of final static entanglement we could achieve for the steady state of 
the static spins. For large coupling strengths, E F (p s ) switches between close to and almost 
1 from round to round depending on the success of splitting. Thus in this case, a threshold 
criterion for p would be that the state of the static pair spends a certain fraction of its time 
in states with E F close to 1, as shown in Fig. |6(b)| In that case, the average E F « 0.88 for 
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11 



7] =95%. In the weak coupling regime, the steady state E F of the static spins is large and 



fluctuates slightly for sufficiently high success rates, as shown in Fig. |6(c)[ The condition is 
that the amount of reduction in static entanglement by a single unsuccessful event should be 
smaller than the accumulated entanglement built up from a certain number N of consecutive 
successful splittings. In this case then r] = jM^, and we find for example, when 6 ~ 0.1, for 
7] ~ 99.9%, then very high entanglement can be maintained (see Fig. |6(c)| ). 
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Figure 6. Plots for the evolutions of Ep for the static spins: (a) of the attractor state pi after 
one round of interactions with the flying singlet pair entering the same lead, against various 
6 E [0, f ]; (b)-(c) examples of the corresponding steady states against number of rounds n 
when no other errors are present. Once the steady state is achieved, the success rates required 
could be lowered. 



In either case, for the current resource to be useful in terms of our scheme, the success 
rates for splittings should improve from those in lfT3l[l4| to at least 90%. Stronger currents 
will also make the scheme more robust to the various decoherence sources. 



6. Other Models 



The above analysis was adapted to a Heisenberg exchange model by replacing the effective 
Hamiltonian in Eq. 1 with 



1 cy \ X W X ' 



v y 



a s Ja f J] 



(17) 



We then obtained similar results for the Heisenberg exchange model; the only difference 
is that this time the singlet state is the only attractor for Eq. [9l and the condition of similar 
coupling strengths is now replaced by the requirement of 6i(k) + 2 (k) ~ n in each round of 
interaction (see Fig.[7]in comparison with Fig. [3]). Given this understanding, the results on the 
decoherence effects as well as the splitting efficiency are also similar. 

The Heisenberg and XY coupling models can be regarded as examples of a more 
general anisotropic coupling E71 . and both can be realized in suitable materials ll27l l28l . 
For example, the Heisenberg exchange model can effectively describe the spin evolutions 
of forward scattering electrons: The evolution of the total wavef unction of the scattering 
electrons is determined by the Schrodinger equation with Coulomb repulsion terms, together 
with the Pauli principle the spin evolution is determined ll28l . 
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Figure 7. Contour plots for the single attractor state of the static spins in the Heisenberg 
exchange model: (a) Fidelity with respect to the singlet state as a function of 6i e [0, it]; (b) 
Ep as a function of Qi £ [0,ir]. 



7. Conclusion 

We have shown that distributed and static spin entanglement can be generated from a source 
of entangled current and weak, passive, Coulomb interactions. The entanglement generated 
is robust to various error sources. We therefore anticipate that spin entangled currents can be 
utilized in the way we have proposed in a wide range of experimental systems. 
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